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ABSTRACT 


The  numerical  solution  of  the  partial  differential  equation  in  the  in¬ 
compressible  turbulent  boundary  layer  has  been  obtained  for  step 


pu,C  Vr/2 
1  P  f 


and  for  Prandtl  numbers  0.7,  1  and  7.  The  Schmidt  method  of 


integration  was  used  and  the  integration  was  carried  out  on  a  Ferranti  Pegasus 
Digital  Computer.  A  method  has  been  developed  to  apply  this  numerical  solution 
for  obtaining  surface  and  fluid  temperature  for  the  case  of  arbitrary  distribution 
of  heat  flux  at  the  surface. 


Simultaneously,  approximate  equations  for  calculating  heat  transfer  over  a 
flat  plate  with  arbitrary  heat  flux  are  derived  for  laminar  and  turbulent  flow. 

To  verify  the  theoretical  solution,  experiments  were  made  in  which  the 
concentration  profiles  of  the  injected  gas  (carbon  dioxide)  at  different  stations 
were  measured  when  the  pipe  had  another  gas  (air)  flowing  turbulently  through  in 
the  axial  direction,  and  the  gas  injected  was  passed  through  a  porous  section  of 
the  tube  wall. 
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1.1.  The  Problem 


The  initiation  of  this  investigation  was  the  idea  that  there  might  be 
important  advantages  in  the  construction  of  a  solid  propellant  grain  by  the 
assembly  of  alternate  "pineapple  rings"  of  different  propellants.  This  problem 
was  therefore  investigated  at  some  length  both  theoretically  and  by  means  of 
an  idealised  experiment. 

In  considering  the  operation  of  such  a  grain,  the  first  problem  is  that  of 
determining  the  manner  of  mixing  the  gases  deriving  from  the  alternate  slices. 

This  problem  idealises  to  the  problem  of  finding  mass  transfer  coefficient 
and  the  concentration  in  the  boundary  layer  for  an  arbitrary  mass  diffusing 
surface.  Heat  transfer  and  mass  transfer  are  similar  processes  with  similar 
differential  equations.  Therefore,  in  the  case  of  heat  transfer,  this  is  a  problem 
of  finding  heat  transfer  and  temperature  profile  in  the  boundary  layer,  for  non- 
isothermal  surfaces.  The  present  research  is  a  contribution  to  the  solution  of 
this  problem. 

In  the  uniform  surface  temperature  case,  analytical  predictions  of  the  local 
heat  transfer,  for  laminar  flow  over  a  flat  plate  and  wedge,  have  met  with 
considerable  success.  When  the  flow  in  the  boundary  layer  is  turbulent,  then 
the  momentum  and  energy  equations  for  the  laminar  boundary  layer  apply  with 
the  addition  of  the  eddy  transport  properties  to  the  viscosity  and  thermal 
conductivity.  However,  these  eddy  transport  properties  are  not  known  fundamentally 
so  that  the  available  solutions  are  of  the  analogy  type,  relating  the  heat  transfer 
to  the  friction. 

Considerable  attention  has  been  directed  to  the  non-uniform  wall  temperature 
case,  with  notable  contributions  by  Lighthill'^,  Leveque^^,  Rubesin^^, 

Eckert^®\  and  Spalding^,  etc.  for  the  laminar  case,  and  by  Seban'*^,  Rubesin'29), 
etc.  for  the  turbulent  case.  Reviews  of  such  solutions  are  given  by  Tribus  and 
Klein^,  Spalding  and  Pun'^,  and  Shah^D. 

Although  there  are  numerous  methods  availble  for  calculating  heat  transfer 
over  the  bodies  of  arbitrary  shape  with  arbitrary  surface  temperature,  there  is 
no  method  available  for  calculating  the  temperature  within  the  boundary  layer. 

1,2.  Present  Work 

In  the  present  investigation,  the  theoretical  approach  to  this  problem  has  been 
made  in  the  field  of  heat  transfer,  and  the  experimental  approach  is  made  in  the 
field  of  mass  transfer. 

On  the  theoretical  side,  the  numerical  solution  of  the  partial  differential 
equation  in  the  incompressible  turbulent  boundary  layer  has  been  obtained  for 
q" 

step  - — —  ■  and  for  Prandtl  numbers  0.7,  1  and  7.  The  Spalding 

P  ui  Cp  Cj/2 

boundary  layer  velocity-law  was  assumed,  and  the  Schmidt  method  of  integration 
used.  Integration  was  done  on  a  Ferranti  Pegasus  Digital  Computer.  A  method 
has  been  given  for  applying  this  solution  to  the  problem  of  obtaining  surface  and 
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fluid  temperatures  for  the  case  of  arbitrary  distribution  of  heat  flux  at  the  wall. 

Simultaneously,  approximate  equations  for  calculating  heat  transfer  over 
a  flat  plate  with  arbitrary  heat  flux  are  derived  for  laminar  and  turbulent  flow. 

These  approximate  equations  are  very  simple  compared  with  similar  equations 
of  Reynolds  et  al^®^,  Rubesin^^  etc.  ,  and  still  predict  fairly  accurate  heat 
transfer  coefficients. 

To  verify  the  theoretical  solution,  an  experimental  rig  was  designed  in  which 
the  concentration  profiles  of  the  injected  gas  (carbon  dioxide)  at  different  sections 
were  measured  when  the  pipe  had  another  gas  (air)  flowing  turbulently  through  in 
the  axial  direction,  and  the  gas  injected  was  passed  through  a  porous  section  of 
the  tube  wall.  The  concentration  of  the  carbon  dioxide  was  measured  by  an  Infra¬ 
red  Gas  Analyser. 

Designing  and  manufacturing  the  experimental  rig  took  a  long  time.  After 
many  modifications,  the  experimental  results  obtained  were  reliable  and  repeatable. 

Experimental  results  showed  higher  diffusion  when  compared  with  the 
theoretical  solution.  In  consequence,  the  experimental  concentration  of  carbon 
dioxide  at  the  surface  was  lower  by  25%  of  the  theoretical  concentration.  It  seems 
probable  that  most  of  the  discrepancy  is  due  to  the  application  of  two-dimensional 
theory  to  the  three-dimensional  case  of  pipe  flow. 

2.  Review  of  the  Previous  Work 


2.1.  Introduction 


In  the  period  of  the  last  ten  years,  research  in  the  field  of  heat  transfer  has 
increased  at  a  tremendous  rate,  and  numerous  investigations  are  directed  to  the 
problem  of  "heat  transfer  from  non -isothermal  surfaces". 

The  method  for  determining  heat  transfer  for  non-isothermal  surfaces  is 
similar  to  the  methods  used  in  determining  the  deflection  of  beams  subjected  to 
arbitrary  load  distributions.  The  energy  equation  of  the  boundary  layer  is  linear 
in  the  fluid  temperature,  if  the  fluid  properties  are  assumed  to  be  constant.  This 
allows  the  super-position  technique  to  be  employed.  Rubesin  has  shown  that  the 
heat  transfer  rate  for  an  arbitrary  temperature  variation  can  be  determined  by 
superimposing  a  number  of  "step  temperature  distributions",  so  that  the  summation 
of  the  steps  yields  the  actual  variable  temperature  distribution,  and  the  heat  transfer 
at  any  point  is  equal  to  the  sum  of  the  heat  transfer  rates  attributed  to  all  steps 
upstream  of  the  point  in  question. 

All  the  approximate  methods,  instead  of  satisfying  the  differential  equation 
for  every  fluid  particle,  satisfy  boundary  conditions  near  the  wall  and  the  region 
of  transition  to  the  external  flow,  together  with  certain  compatibility  conditions. 

In  the  remaining  region  of  the  fluid  in  the  boundary  layer,  only  an  average  value 
over  the  differential  equation  is  satisfied,  the  average  being  taken  over  the  whole 
thickness  of  the  boundary  layer.  Such  mean  values  can  be  obtained  from  the 
momentum  theory  and  the  law  of  the  conservation  of  energy.  These  are,  in  turn, 
derived  from  the  equations  of  motion  and  energy  by  integrating  over  the  boundary 
layer  thickness.  These  integral  equations  of  motion  and  heat  flux,  with  assumed 
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velocity  and  temperature  profiles  in  the  boundary  layer,  form  the  basis  of 
the  approximate  methods. 


In  all  cases,  if  the  wall  temperature  is  given, 
from 


qw(x) 


t 

/ 


?=x 

?=0 


l(C,x)  d  Tw(0 


the  heat  flux  is  to  be  calculated 


(2.1) 


Here  the  Kernel  function  h(f,x)  is  the  heat  transfer  rate  at  position  x  due  to  step 
temperature  rise  at  the  position  If  the  heat  flux  is  given,  the  wall 

temperature  is  to  be  calculated  from 

r 

°-2' 

Here  the  Kernel  function  g,  f  .is  the  wall  temperature  at  x  due  to  unit  heat 
flux  at  5.  '5'X, 


It  should  be  noted  that  integral  of  equations  2. 1  and  2. 2  muBt  be  taken  in  the 
"Stieltjefl"  sense  rather  than  in  the  ordinary  "Riemann"  or  "Area"  sense. 
Various  investigators  have  obtained  analytical  expressions  for  the  "Kernel" 
functions,  and  all  these  expressions  can  be  divided  into  four  main  groups 
depending  on  their  applicability. 


2.2.  Uniform  Main  Stream  Velocity  -  Laminar  Flow 


(23) 

Leveque  assumed  a  linear  velocity  profile  independent  of  x,  and  solved 
the  differential  energy  equation.  He  obtained 
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<C,x) 


k(Pr)3 

3(1)]“ 


-i 

(x  -  0 


(2.3) 


(15) 

Rubesin  assumed  linear  velocity  and  temperature  profiles  and  solved 
the  integral  energy  equation  by  further  assuming  that  the  thermal  boundary 
layer  varies  proportionally  to  the  momentum  thickness. 


He  obtained 
h 

.(9) 


0.304k 


(C.x) 


1  1  r-  w  3  ^ 

Pr’  [*-$•] 


‘I 


(2.4) 


Eckert"'  assumed  cubic  velocity  and  temperature  profiles  and  solved  the 
integral  energy  equation.  He  obtained 


(f.x) 


0.33k  _  £  _ 

-  Pr  Re 

x  x 


L  j—  f  .3.  n 

2  i  -  (V 

X  x 


(24) 


Klein  and  Tribus'"  "■  have  shown  a  mathematical  manipulation  procedure 
by  which  the  value  of  ^  can  be  obtained  from  the  known  function  h^  ^ 


Thus  the  values  of  the  function  g.„  .  obtained  from  the  equations  of 

,  (23)  „  ,  .  (15)  ,  ,  ’%) 

Leveque  ,  Rubesin  and  Eckert'  '  are 
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g(C,x> 

-l  .4  _ 

2  Pr  3  /  p  y  /du\ 

3k  (1)1  V W  [Wy=0. 

(x  -  C)  5 

(2.6) 

g(C,x) 

1  Pr  n  i 

8(i)l  (|)1  0.304k  x 

(2.7) 

and 

g(C,x) 

1  Pr"*  „  -£  r. 

6(J)1  (|)1  0.33k  x 

**  r 

(2.8) 

respectively. 

The  present  authors  have  assumed  cubic  velocity  and  temperature  profiles 
and  have  obtained  a  very  simple  equation  for  an  arbitrary  heat  flux  problem  by 
directly  integrating  the  energy  equation  for  a  step  heat  flux  case.  Detailed 
derivation  is  given  in  Chapter  3. 


2.3.  Uniform  Main  Stream  Velocity  -  Turbulent  Flow 

Rubesin^^  assumed  velocity  and  temperature  profiles  as  fo 
power  law,  and  solved  the  integral  energy  equation.  He  obtained 


0.0288k  „  i  ,0.8 
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and 
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.  .  .  (2.10) 


Seban  followed  the  same  procedure  as  Rubesin,  except  that  he  assumed 
a  linear  temperature  profile  in  the  laminar  sublayer.  He  obtained 


0.0289k  „  1/9  „  0.8 

a,.  .  =  -  Pr  Re 

(C.x)  x  x 


1  -  A 
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g(C,x)  =  (8/9)'.  <1/9)1  0.0289k 


,(14) 
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Maisel  and  Sherwood  carried  out  experimental  measurements  on  a  mass 
transfer  apparatus  and  obtained  an  empirical  equation  "best  fit"  to  their  data. 
For  air:- 


0.035k 


(C.x) 

,(10) 
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(2.13) 


Reynolds  et  al  solved  the  integral  energy  equation  by  assuming  velocity 
and  temperature  profiles  as  following  1/7  power  law  and  turbulent  Prandtl  number 
equal  to  unity.  They  obtained 
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For  the  step  heat  flux  problem,  the  present  authors  have  derived  a  simple 
analytical  expression  by  assuming  a  seventh  power  law  for  velocity  and  temperature 
profiles  and  directly  integrating  the  energy  equation.  The  detailed  analysis  is 
given  in  Section  3. 

2.4.  Arbitrary  Main  Stream  Velocity  -  Laminar  Flow 
(5) 

Lighthill  solved  the  energy  equation  in  Von  Mises  form,  assuming  a  linear 
velocity  profile  in  the  thermal  boundary  layer.  He  obtained 
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(2.17) 


(19) 

Bond  also  assumed  a  linear  velocity  profile  and  solved  the  differential 
equation  for  wedge  flow.  For  the  wedge  flow  (ut  =  Axm)  he  obtained 
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is  a  dimensionless  velocity  gradient 


(2.21) 


(7) 


Ambrok'  '  has  derived  an  approximate  method  for  calculating  the  heat  transfer 
coefficient  by  assuming  that  a  relation  of  the  type  Nu  =  A  Re  (Ag)n,  which  is  valid 
for  a  flat  plate  with  a  uniform  surface  temperature,  and  is  also  valid  for  a  problem 
where  main  stream  velocity  and  surface  temperature  vary  arbitrarily.  With  this 
assumption  he  has  solved  an  integral  energy  equation.  He  obtained 
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This  equation  is  very  simple  but  it  is  found  to  give  low  values  of  heat  transfer 
when  compared  with  the  exact  solution  of  Eckert  for  wedge  flows. 

(8) 

Spalding  has  improved  on  Lighthill's  method  by  a  correction  which  accounts 
for  the  departure  from  linearity  of  the  velocity  profile  within  the  thermal  boundary 
layer,  and  which  comprehends  the  influence  of  Prandtl  number,  pressure  gradient, 
body  forces  and  non-coincident  start  of  velocity  and  thermal  layers. 

In  this  method,  the  momentum  thickness  6  ,  of  the  velocity  boundary  layer 
is  to  be  evaluated  first  by  a  procedure  similar  to  those  of  Walz^6)  an(j  Xhwaites^5). 
With  6  (x)  evaluated  6  (x)  is  to  be  obtained  from  a  table  or  graph  of  6/6  vs. 

6*  / du  ‘ 

—&■  (  — !  )  ■  Finally,  the  heat  transfer  coefficient  to  be  obtained  from 
»  \dx/ 

l  r  -  i  -l 

/  6  ,2  r  ,  u  *  ?  /u  2 

=  [£)  [6'41*  j  [ r-J  dx+aj  Fdx  (2.23) 

a  a 

6  b  du 

where  F  is  a  graphically  presented  function  of  the  argument  -  -  . 

v  dx 

(4) 

Schuh  selected  a  dependent  variable,  which  is  a  function  of  the  ratio  of 
heat  flow  across  the  whole  boundary  layer  to  the  temperature  gradient  at  the  wall 
in  a  suitable  dimensionless  form,  and  then  integrated  the  momentum  and  the  energy 
equations. 

Schuh's  method  involves  similar  steps  to  those  of  Spalding's  method,  namely: - 

1.  The  determination  of  62  by  the  Walz-Thwaites  technique. 

6 2  du 

2.  fi  and  Z,  are  the  functions  of  -z-  — r-5-  and  to  be  read  from  the  graph. 

V  dx  b  r 

3.  Calculation  of  n  from  equation  2. 24 


tJ  F 


(2.24) 


4.  Calculation  of  functions  G  and  P  from  equations  2.25  and  2.27  by  an 
iterative  procedure. 

G  =  0.57  (/3  +  0.205)0'104  A°'37+0-06^  (2-/3)2  P1_3r  (2.25) 


whe  re 


A  =  1  +  <2  -  0)  n 


(2.26) 


and  values  of  r  are  to  be  read  from  his  table. 


rjJ  a1'6  F 

M  [ 

7  [TLl 1  ' 

a  [_(2-/3)x 


(T  -  T  )  G(2  -  $) 

W  1 


T  [u,x] 


i  g3/2  ('i 


(T  -  T  )" 

w  ' 


(2.27) 
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5.  Heat  transfer  coefficient  then  to  be  calculated  from 


v 


(2.28) 


Seban  and  Drake  solved  the  momentum  equation  by  a  procedure  similar 
to  that  of  Eckert,  and  obtained  heat  transfer  rate  by  assuming  tha.t  it  is  related  to 
the  momentum  thickness  in  the  same  manner  as  on  the  wedge. 


This  method  involves  similar  steps  of  those  of  Spalding  and  S»chuh. 

6*  du 

1.  The  determination  of  &a  and  the  auxiliary  function 

2.  Evaluation  of  functions  z(  and  (1  from  his  graph. 


3. 


The  rate  of  heat  transfer  at  the  wall  for  a  particular  case,  wlien  wall 
temperature  is  a  power  function  of  length,  e.g.  (T  -  T^  ■  -A  xn,  is  given  by 


.  *  ,  .  n  *  * , 

q  =  -  k  A  x  g.  .  — 
Mw  *(o)  6 

2 


(2.28) 


or 


(2.30) 


where  g^  is  a  temperature  gradient  at  the  wall,  g^  is  a  function  of  0,  n  and  Pr. 

Seban  has  suggested  the  exact  solutions  of  Eckert,  Levy  etc.  to  b«  used  to  obtain 
the  value  of  gj^  for  different  /},  n  and  Pr. 

By  superimposing,  the  heat  transfer  coefficient  can  be  deterimined  for  a 
variation  of  the  surface  temperature  of  the  type 


n,  n  j 

(T  -T)  *  Ax  +  A  x  +  Ax 
w  <  i  r  i 


(2.31) 


2.5.  Arbitrary  Main  Stream  Velocity  -  Turbulent  Flow 

(16) 

Very  recently  Spalding'  '  has  suggested  a  single  analytical  formula  for 
velocity  distribution  throughout  the  whole  of  a  "Universal"  turbulent  boundary 
layer  in  the  form  of  y+(u+)  instead  of  the  usual  form  u+(y+).  Wi-fch  this  form  of 
expression,  it  is  possible  to  obtain  numerical  solutions  of  the  par-tial  differential 
equation  for  arbitrary  surface  temperature  case  and  for  the  arbitrary  heat  flux 
case.  These  solutions  cover  both  the  laminar  and  turbulent  flows  . 


( 1 6) 

Spalding'  ;  has  solved  the  integral  form  of  the  energy  equation  for  step 
surface  temperature  and  for  unity  Prandtl  number.  He  assumed  bl  linear 
temperature  profile  of  the  form 


u+  <  6  ;  6 
u+  >  6  ;  0 


6 


3 


0 


(2.32) 
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(17) 

Muralidharan  has  solved  the  differential  equation  for  step  surface 
temperature  and  for  Pr  -  0.7,  1  and  7.  He  used  an  analytical  formula  of 
Deissler^^  for  velocity  distribution  in  the  boundary  layer  and  obtained  the 
numerical  solution  on  an  analogue  computer. 

( 1 8) 

Kestin  and  Persen  used  the  digital  computer  to  obtain  the  numerical 
solution  of  the  partial  differential  equation  for  step  surface  temperature  and 
for  Pr  =  1 .  They  used  Spalding's  velocity  law. 

The  authors  have  solved  the  partial  differential  equation  for  arbitrary  heat 
flux  problem  and  for  Pr  =  0.7,  1  and  7.  The  detail  of  the  numerical  solution  is 
presented  in  Section  4. 

3.  Heat  Transfer  over  a  Flat  Plate  with  Arbitrary  Heat  Flux.  Approximate 
Analysis 

3.1.  Laminar  Flow 


m(x-a) 


(  a  ) 

•  J 

. 

<  - ? 

*  A 

Fig.  3.2 


For  a  step  heat  flux  (Fig.  3.1),  total  energy  at  any  section  x  is:- 

qMx  -a)  =  pUi  Cp(Tw  -T)6  J  <1  -  9 )  d  (f)  (3.1) 

o 

For  assumed  cubic  velocity  and  temperature  profiles 


Equation  3.1  becomes 
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q  (x  -  a)  =  p  u  C  {T  -  T 
^w  i  p  w  ‘ 


i’  6  (j)  '  ik  (s)  ]  (3>4) 


We  are  treating  the  case,  that  the  thermal  boundary  layer  is  thinner  than 
the  velocity  boundary  layer,  that  is  (A/ 6)  <  1,  The  second  term  in  the  right- 

hand  expression  of  equation  3.4  is  then  small  compared  with  the  first,  and  can 
be  neglected.  The  heat  transfer  at  the  wall  is  given  by 


h(T  -  T  )  =  -  k 
w  * 


ill  _  1 

_Hy=0  =  2 


k(T 


T) 


(3.5) 


For  a  laminar  boundary  layer,  velocity  boundary  layer  thickness  is  given  by: 

4.64  x 


v'r! 


(3.6) 


Thus  for  a  step  heat  flux,  from  equations  3.4,  3.5  and  3.6 


(T  -  T  )  = 
w  1 


1  i 

2.395  b(Re)2  (Pr)3 

p  u  C 
>  P 


[>  -fc 


and 


Nu 


/Re 


X  =  0.418  (Pr)^ 


l*  -  wr 


(3.7) 


(3.8) 


Where  b  is  the  magnitude  of  the  step  heat  flux. 

Similar  equations  of  Eckert^  and  Rubesin^15^  for  a  step  heat  flux  are 

_T)  ,  2.2  b  (H.)*  IPD1  ,  „ 
w  t  p  u  C  r 


and 


,  ,  2.405  b  (He)*  IPrl*  ,  „  4;„ 
w  i  pu  C  r 


(3.9) 


(3.10) 


respectively.  It  can  be  seen  that  equation  3.7  is  relatively  simple,  for  it  does 
not  contain  the  incomplete  beta  function.  For  any  arbitrary  heat  flux,  the  wall 
temperature  and  the  heat  transfer  coefficient  can  be  calculated  by  superimposing 
small  steps.  Thus  for  any  arbitrary  heat  flux 


(  T  -  T  )  = 


1  2 
2.395  ReT  Pr* 

p  u  C 
1  P 


S=x 


r=o 


(1  -  S/x) 


i  d  ^(g) 

dC 


df 


(3.11) 


and 


Nu 
_ x 


0.418  (Pr)3  q* 


(x) 


/, 


?=0 


1  d  ^w(f) 

(1  -  C/x>*  dC 


(3.12) 
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It  should  be  noted  that  the  integral  of  3.11  and  3.12  must  be  taken  in  the 
"Stieltjes"  sense.  For  instance,  for  a  delayed  ramp  heat  flux  (Fig.  3.2): 

i  2  4/'j 

t  /n.lJ  ~ t  ’  '  ° 


(T  -  T  )  = 
w  1 


1 .797  (Re)2  (Pr)3  mx 
p  u  C 


[  1  -  (a/x) 


(3.13) 


=  0.557  (Pr)5  1  -  (a/x) 


(3.14) 


3.2.  Turbulent  Boundary  Layer 

For  a  step  heat  flux,  in  the  turbulent  boundary  layer,  we  assume  velocity 
and  temperature  profiles 


r  i=  (y/6) 


(3.15) 


T  -  T  . 

w  ,  , .  v  1  /  7 
T  .  T  =  (y/A) 


(3.16) 


Then  equation  3.1  becomes 


q*  (x  -  a)  -  p 
w 


u  C  (T  -  T  )6  4 

i  p  w  i  72 


*  & 


(3. 17) 


The  heat  transfer  momentum  analogy  derived  by  Reynolds  et  al  ,  for  step 
heat  flux  is:  - 

cr  /»  :1/7 


(T  -  T  )  pu  C 
w  i  ip 


T 


where  the  Prandtl  number  dependence  factor  is  applied  afterwards. 
For  a  turbulent  boundary  layer 


(3.18) 


4-  =  0.0296  (Re)'°‘2 


(3.19) 


-  =  0.37  (Re)'0,2 

x 


(3.20) 


Thus  for  a  step  heat  flux,  from  equations  3.17,  3.18,  3.19  and  3.20  :  - 


32.95  b  .0.2  .0.4.,  .  ,1/9 

(T  -  T.)  =  —7: -  (Re)  (Pr)  (1  -  a/x) 

w  1  p  C  u 
P  » 


(3.21) 


0.03035  (Pr)0,8  (1  -  a/x)'1^9 


(3.22) 


where  b  is  the  magnitude  of  the  step  heat  flux.  The  similar  equation  of  Renolds  et  al 


for  a  step  heat  flux  is:  - 
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(T  -  T  )  =  32.42 
w  i 


p  C  u 
P  ’ 


(Re)0'2  (Pr)0'4  I  (1/9,  10/9) 


(3.23) 


Equation  3.21  is  relatively  simple,  for  it  does  not  contain  the  incomplete  beta 
function.  For  any  arbitrary  heat  flux,  the  wall  temperature  and  the  heat 
transfer  coefficient  can  be  calculated  by  superimposing  small  steps.  Thus  for 
any  arbitrary  heat  flux:- 


(T  _  T  )  _  32.95  (Re)°  2  (Pr)°'4 

w  i  p  C  u 


and 


S=x 

I  (l  -  C/x) 
JC=0 


1/9  dqw(g) 


dC 


d? 

(3.24) 


Nu 


0.03035  (Pr) 


OR 


qw(x) 


(Re) 


0.8 


C--x 


f-o 


(3.25) 


0  -  C/x)1/9  d-Sl>dC 


In  equations  3.24  and  3.25  the  integral  must  be  taken  in  the  "Stieltjes"  sense. 
Taking  a  delayed  ramp  heat  input  (Fig.  3.2)  as  an  example:- 


(T  -T,l  ■  29-65  mJ,lP,),'  V|0f  ,,..,,,10 /• 


and 


Ku 


(Re) 


0.8 


P  C  u 
P  1 


0.03373  (Pr)0-6  (1  -  a/xf 1/9 


(3.26) 


(3.27) 


Results  of  the  above  equations  are  compared  with  Reynolds  et  al 
Eckert'®',  and  Rubesin^^,  for  both  laminar  and  turbulent  boundary  layers, 
and  are  plotted  in  Figs.  3.3  and  3.4. 


For  the  laminar  case,  the  results  of  the  above  equations  agree  very  well 
with  that  of  Rubesin'*  ,  and  not  with  Eckert^®*.  The  disagreement  with  Eckert 
is  merely  in  a  numerical  factor.  Eckert's  (T  -  T,)  values  are  8%  lower  than  those 
of  equation  3.7. 


fi'or  the  turbulent  case,  the  results  agree  within  1.5%  to  that  of  Reynolds 

l0>. 
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3.3.  Approximate  Analysis  assuming  Linear  Temperature  Profih 
For  a  step  heat  flux,  total  energy  at  any  section  x  is: 


q "  (x  -  a)  = 
w 


PC  u  (T  -  T, )  dy 
P 


St  (x  -  a)  u  = 


r 


6  dy 


With 


T  w  _  (x  -  a) 
p  v  M  2 


c-t 

it  u  y  1C 

X  L*  „  _il_  I- 

v  p  v  N 

_u_  _  u  _1 

I  rw  Ut  I  Cf 

N  p  V  ”5" 


and 


4X1  =  e+ 

du+ 

equation  3,29  simplifies  to: 

6 


St 


f 


I  u+  6  e+ 


du 


Assuming  a  linear  temperature  profile 

ill 

6 


e  =  l 

and  knowing  that 


•\ 


r  ■-(“)♦„ 


Equation  3.34  simplifies  to: 


x 

6 


6  /  .  +\ 

I  u+  (1  *  S’ )  e+  du+ 


(16) 

From  Spalding's'  '  velocity  profile:- 


e+  =  1  + 


k'  r  k'u+ 

E  |_e 


1  -  k  'u 


»„+ 


with 

and 


k'  =  0.407 

=  0.0991 
E 


(k'u+)a 

2'. 


(k'u*)3  ] 

3-.  J 


(Pr  =  1) 

(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 
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From  equations  3.37  and  3.38 


_  =  o.  +  _ 

6  6  U  e'  E 


r  k  6 
e 


a  -  4- )  +  -„u  +  4-> 

k  &  k s  k  6 


It'S5  f1  *  3  k'Vj 

‘  IT  LZ  +  20  +  IF  J 


(3.39) 


Equation  3.39  is  compared  with  previous  approximate  analyses  for  laminar  and 
turbulent  flows  obtained  by  assuming  cubic  and  seventh  power  profiles  repsectively. 
(Eqs.  3.42  and  3.45).  Comparison  is  shown  in  Fig.  3.5  and  agreement  is  found 
to  be  very  satisfactory. 


For  Pr  =  1 ,  equation  3.  8  of  laminar  flow  is: 
—  =  0.418  (1  -  9/xf* 


vrRe 

Knowing  C 


0.332 

fRe 


for  laminar  flow,  equation  3.40  simplifies  to: 


J 


I~  =  0.603  x+'  ^ 

Cf 

~T 


For  Pr  =  1 ,  equation  3.  22  of  turbulent  flow  is: 


0.8 


=  0.03035  (1  -  9 /x) 


with 


Nu_ 

Re 

Cf  -0  2 

~  =  0.0296  Re 


-1/9 


equation  3,43  simplifies  to: 

=  0.1452  x+ 

C, 


,-1/9 


N  2 


(3.40) 

(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45) 
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4.  Numerical  Solution  for  Arbitrary  Heat  Flux 


4.1.  Theory 

( 1 6) 

By  use  of  the  Von  Mises  transformation,  Spalding'  has  reduced  the  energy 
equation  for  a  turbulent  incompressible  boundary  layer 


(this  implies  that  shear  stress  is  independent  of  y) 


Pr 


Further,  using  a  new  variable  £  defined  by  d&  =  — —  du4 

«+ 

equation  4,2  can  be  reduced  to 


3T  _  _1_  d1  T 

dx+  "  u+a+  ae* 

with  boundary  conditions 


T  =  0  at  £  =  «>  ) 

and  T  =  0  at  x+<  0  ) 


(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 


In  this  section  the  temperatures  T  and  are  relative  to  the  fluid  temperature 
An  additional  boundary  condition  in  the  present  solution  is 


Constant 


(4.11) 
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(16) 

Spalding'  has  derived  a  single  expression  for  the  distribution  of  the 
velocity  in  the  universal  turbulent  boundary  layer  in  the  form  of  y+  (u+) 
instead  of  the  usual  form  u+(y+).  It  is:  - 


where  k* 
and  1 / E 


i  j  kV 

i[e 


0.407 

0,0991 


1  -  kV 


(kV)1  0cV)3  (kVr  i 

2’.  "  3*.  "  4>.  J 


(4.12) 


From  equation  4. 12  an  expression  for  e+  can  be  obtained,  using  equation  4.6. 
This  gives 


II 


k  u+  i  + 

e  -  1  -  k  u 


(kV)3 

3! 


(4.13) 


A  further  assumption  that  the  turbulent  Prandtl  number  is  unity  (e^  =  e)  leads 
to  an  expression  for  a+ 


+  1  k'  j  k'  u^  ,  ,  /  + 

'  =  Pr  +  E  L 6  -1  -ku  - 


(k'u+)a 


2! 


(k  uf) 
3'- 


(4.14) 


Therefore  the  quantity  (u+  <*+),  for  various  Prandtl  numbers,  can  be  obtained  as 
a  function  of  4  .  This  reduceB  equation  4.9  to 


3T  1  af_T 

9x+  *  Ad  ae* 


(4.15) 


The  values  of  f(4)  for  Pr  =  0.7,  1  and  7  are  shown  in  Fig.  4.6. 

For  the  uniform  surface  temperature  case,  various  investigators  have 
solved  equation  4. 15.  Spalding*1”'  has  solved  it  approximately  for  Pr  =  1 ,  by 
use  of  the  energy  integral  equation  and  assuming  a  temperature  profile  in  the 
boundary  layer.  Muralidharan'1^  has  solved  it  on  an  analogue  computer  for 
Pr  =  0.7,  1  and  7.  Kestin  and  Persen^1®^  have  solved  it  on  a  digital  computer 
for  Pr  =  1 . 

In  the  present  paper,  a  solution  of  equation  4. 15  has  been  made  for  uniform 
at  t*le  wa^  anc*  ^or  1>r  =  1  arKl  ^  ■ 

4.2.  Numerical  Solution 


It  will  be  seen  that  equation  4.15  is  very  similar  to  the  heat  conduction 
equation,  and  may  be  solved  by  the  finite  difference  method  of  E.  Schmidt. 


The  equation  4.15  can  be  written  in  finite -difference  form  as 


T(x+  + Ax+,4)’  V,  4) 


I  T 

f(£)  (<£)*  [_  <x*  S 


+  T(x+,  4  -A4) 


■V.n] 

(4.18) 
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From  equation  4. 16,  T^+  +  ^x+)  can  be  calculated  from  which 

in  turn  can  be  calculated  from  T^x+  _  and  so  on.  The  grid  for  the  finite 

difference  scheme  is  shown  in  Fig.  4.7. 


From  Fig.  4.  7  it  can  be  seen  that  all  the  values  of  temperature  can  be 
calculated  for  T^+  +  ^x+j  from  Tjx+j  except  two  end  values,  one  at  the 

wall  and  the  second  at  t  =  ^max‘  *a*ter  difficulty  was  met  by  selecting 

4  so  high  that  T  there  was  zero.  The  value  of  T  at  the  wall  was  obtained 

max 

by  adding  the  wall  temperature  gradient  x  A£  to  the  value  of  the  temperature 
at  (x+  +  Ax+,  A 4). 


To  obtain  the  numerical  results  in  solving  equation  4.15,  the  complete 


calculations  were  made  for  an  arbitrary  value 


(  8T\ 

\34/4= 0 


=  -  40, 


and  later 


the  results  were  made  dimensionless.  The  number  -40  has  no  particular 
significance. 


= 


all  values  of  4  at  x+  -  0  were  zero.  This  difficulty  was  overcome  by  starting 
the  solution  at  a  very  low  value  of  x+  called  x+  . 


The  starting  value  of  x|  for 


MJr_ 


-  -40  was  selected  such  that  the 


4=0 


temperature  at  all  values  of  4  was  zero  except  at  the  surface.  This  first - 
assumed  value  of  the  surface  temperature  T  was  obtained  from  the  approximate 

(35)  Wl 

analysis  of  Smith  and  Shah  The  final  results  when  checked  by  integrating 

the  heat  flux  showed  that  they  were  not  affected  by  any  error  in  this  assumption. 
The  starting  values  for  various  Prandtl  numbers  are  tabulated  in  Table  4  ,1. 


Pr 

T 

W  1 

x+ 

1 

A4 

0.7 

19.30 

0.050 

0.5 

1.0 

19.41 

0.025 

0.5 

7.0 

38.50 

0.004 

1.0 

Table  4.1 


The  solution  by  the  method  of  E.  Schmidt  is  stable  only  if  the  condition 


f(4 )  (A  4  ) ' 

is  fulfilled  throughout.  Since  it  was  decided  to  solve  the  equation  up  to 
x+  =  10®,  the  interval  A  x+  was  increased  gradually  as  shown  in  Tables  4.2, 
4.3  and  4.4.  With  the  increase  in  the  interval  Ax+,  it  became  necessary  to 
increase  A4  near  the  surface  (low  values  of  4 )  to  fulfil  the  above  condition. 
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As  the  values  of  f(C)  increase  with  C  the  interval  AC  was  gradually  reduced  for 
the  calculations  away  from  the  surface.  Near  the  surface,  the  temperature 
profile  gradually  became  linear  and  therefore  this  increase  in  AC  near  the  wall, 
with  the  increase  in  the  interval  Ax+  did  not  impair  the  accuracy  of  the  final 
results. 


Interval  x+ 

Ax+ 

Min.  AC 

Max.  AC 

0  to  1 

0.05 

0.5 

0.5 

1  to  10 

0.1 

0.5 

0.5 

10  to  10a 

1 

0.5 

1 

10*  to  10J 

1 

0.5 

1 

104  to  10s 

50 

0.5 

4 

10s  to  104 

500 

0.5 

8 

Table  4.2.  Steps  chosen  for  Pr  =  0,7 


Interval  x+ 

Ax+ 

Min.  AC 

Max.  AC 

0  to  1 

0.025 

0.5 

0.5 

1  to  10 

0.5 

1 

1 

10  to  10* 

1 

1 

1 

lO*  to  103 

4 

1 

2 

103  to  10‘ 

20 

1 

4 

104  to  10* 

25 

1 

4 

10*  to  10‘ 

250 

1 

8 

Table  4.3.  Steps  chosen  for  Pr  =  1 


Interval  x+ 

Ax+ 

Min.  AC 

Max.  AC 

0  to  1 

0.005 

1 

1 

1  to  10 

0.05 

1 

2 

10  to  10* 

0.5 

1 

4 

10*  to  10* 

1 

2 

6 

103  to 104 

5 

2 

8 

104to  10’ 

25 

2 

14 

Table  4.4.  Steps  chosen  for  Pr  =  7 
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,  Form  of  Presentation  in  the  Results 

The  numbers  resulting  from  the  computations  are  temperature  differences 


T  in  ? ,  x+  co-ordinates,  with 


Heat  transfer  coefficients  may 


'9^  =  0 


be  computed  conveniently  from  such  results  in  terms  of  the  "Spalding  function" 
St/  /Cj./ 2  ,  by  the  relation,  deduced  from  equations  4. 4 ,  4.8  and  4.12. 


■  '(If A  „  /  % 


(4.18) 


Thus  St/V~Cj./2  (x+)  can  be  presented. 
Note,  further,  that 


pCu/CJ2  T 
p  ^  f  v 


(4.19) 


Assuming  that  C^./ 2  (x)  and  1  i(x)  are  known  for  a  specific  example, 
x+(x)  maybe  computed  from  equation  4.3.  Hence  from  St/^C./2  (x+),  St(x) 
maybe  computed.  Then  knowing  the  constant  value  of  q^  /  pu^C^C^I  2  ,  Tw 
may  be  found.  Treatment  for  arbitrary  is  given  in  a  later  section. 

For  the  temperature  within  the  boundary  lajer,  it  will  then  suffice  if 
6  (?,  x+)  be  presented. 

A  practicable  presentation  is  thus  of  St/1^C^/2  (x+)  and  0(£,  x+). 

4.4.  Check  by  Heat  Balance 

An  integral  form  of  the  equation  4.15  may  be  derived  by  integrating  it  along 
the  §  co-ordinate.  Thus 


Constant  and  therefore 


d  ”  ...  +  +  _  /  8T 

S*  /  '  “  “  ds  '  -\n).a 

o 

In  our  solution  (  )  =  Constant  and  therefor 

'0?</f=O 

'  -  (*),„  ** 

Jo 

Therefore  from  equations  4.18  and  4.21  we  have 


(4.20) 


J  6 u+  a+ 


df  = 


(4.21) 


(4.22) 
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The  temperature  profiles  obtained  in  the  numerical  solution  were  integrated 
at  various  sections  and  the  integrated  heat  fluxes  were  compared  with  the  Spalding 
function  and  x+  as  shown  in  equation  4.22.  Comparisons  for  various  Prandtl 
numbers  are  tabulated  in  tables  4.5,  4.6  and  4.  7.  The  results  were  found 
reasonably  satisfactory. 


4.5.  Dependence  of  the  Spalding  function  and  8  on  the  Prandtl  number 

The  results  of  the  calculations  in  terms  of  the  temperature  profiles  are 
shown  in  Figs.  4.1,  4.2,  and  4.3.  The  values  of  the  Spalding  function  for  various 
values  of  x+  and  for  Prandtl  number  0.7,  1  and  7,  are  tabulated  in  Tables  4.8, 

4.9  and  4.10  and  plotted  in  Fig.  4.4. 


Having  performed  these  calculations  and  checks,  an  effort  was  made  to  find 
the  dependence  on  Prandtl  number  of  the  Spalding  function  and  the  temperature 
profile  in  the  form 


/  21' 


2  Pr=l 


Pr 


(4.23) 


Results  for  Pr  »  0.7,  1  and  7  showed  that  in  the  laminar  region  (up  to  x+  *  1000) 
the  value  of  n  remains  constant  and  is  equal  to  -2/3,  but  for  higher  x+,  the  value 
of  n  when  obtained  from  the  results  for  Pr  *  0.7  and  1,  changes  from  -2/3  at 
x+  *  1000  to  -0.4  at  x+  *  10®  and  -0.325  at  x+  =  10”,  whereas  from  the  results 
for  Pr  *  7,  the  value  of  n  decreases  from  -2/3  at  x+  =  1000  to  -0.516  at  x+  =  10®. 


In  the  laminar  region  (x+  <  1000)  the  temperature  profile  for  any  Prandtl 
number  can  be  obtained  from 


where 


®(C)  "  V)(Pr=l) 


{ 


f 'x  Pr 


2/3 


(4.  24) 
(4.25) 


Temperature  profiles  for  various  Prandtl  numbers  and  for  x+  *  100  are 
shown  in  Fig.  4.5. 

4.6.  Application  of  the  Solutions  to  the  Case  of  Arbitrary  Heat  Flux  at  the  Wall 


The  working  equation  by  which  the  numerical  solutions  given  in  this  paper 
may  be  applied  is  equation  4.33  below.  A  full  statement  of  its  development  is 
given  in  Sections  (a)  and  (b)  following. 


(a)  Determination  of  temperature  at  the  wall. 


The  numerical  solutions  given  in  the  previous  sections  have  the  boundary 


condition 


( 

\  85. 


=  constant,  x+  >  0.  From  the  definitions  of  the  variables 


«-0 


in  equations  4.8,  4.7,  4.6  and  earlier,  together  with  the  velocity-law  equation  4.12, 
this  boundary  condition  is  the  same  as  q^/  P  u(  C  V  Cf/2  *  constant. 
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That  is,  the  heat  flux  parameter  is  constant. 

This  boundary  condition  is  highly  particular:  usually  problems  will  be  met 
in  the  form  of  a  specified  distribution  of  q*  ,  with  T  (x)  the  quantity  desired 
to  be  known. 


However,  the  treatment  of  the  problem  of  arbitrary  q "  will  be  most  easily 


perceived  after  a  recapitulation  of  the  method  of  determination  of  Tw(x)  for 
constant  q^/  pu,  C  ^Ctl2  . 


Fig.  4 . 4  gives  values  of  the  Spalding  function  St/‘^C^/2  (Pr,  x+).  If 


St/vrcj72 


be  known  at  a  given  x,  then  T^  is  given  by 

Tw  •  s^-2 « (,7$hvp  ) 


(4.26) 


In  Fig.  4.4,  heat  injection  starts  at  x  =  a.  We  set  the  problem,  then,  of 
calculation  of  the  temperature  T^^,  at  x  =  b,  consequent  on  constant 

q*  jp  u  C  '^C./2  between  x  =  a  and  x  -  b. 

'w  '  p  f 

{a,b 


P  VT 


w/P 


a,  b 


dx 


(4.  27 


t+  , 
a,  b 


Then  xT  ,  is  easily  calculable,  knowing  r 


p  and  v  .  From  Fig.  4.4 
ThiB  quantity  is  the  value  of 


(St/^C,/2)  ,  may  be  read  off,  for  x+  -  x+  ,  . 

t  a, D  a, b 

St/  C./2  at  x  =  b  consequent  on  constant  q"  /  p  u  C  J~Cr/2  from  x  *  a  to 

I  w  “  JJ  I 

x  *  b.  Then  at  x  =  b,  the  difference  T^^)  between  wall  temperature  and 
mainstream  is  given  by 


T 


■(b) 


(- 

\  P  u. 


w 

C  VC./2 
P  I 


(4.28) 


It  must  be  particularly  noted  that  equation  4. 28  holds  for  constant 

q *  /  pu1Cp''^Cp2  between  x  =  a  and  x  =  b.  However,  in  the  general  problem, 

q"/pu.C  *cT/2  will  be  a  function  of  x.  Knowing  q"  ,  ,,  u„  ,  and  C„  .the 
Mw'  ’  p  f  _  6  Mw(x)’  «(x)  f(x) 

dependence  of  q^/  putC^  ^ C^l 2  on  x  may  be  computed.  The  temperature  at 

the  wall  at  x  =  b  is  then  given  by  the  Stieltjes  integral. 

b 


'  w(b) 


’  / 


1 


o  (St/Scp) 


a,b 


(p“, CpVc^Tl) 


(4.29) 
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Equation  4.29  is  the  working  equation  for  the  determination  of  wall  temperature. 

In  equation  4.29,  (St/^C^/2)a  ^  is  the  value  of  St/^Cf/2  at  b.  consequent  on 

•  h  f — “ — ■  *  /  Q  vv  \ 

constant  q  P  u,C  vCr/2  between  x  =  a  and  x  =  b.  d  - — — 3= -  )  is  the 

-  W  ’  p  f  \PUC  /cf/2 / 

increment  of  heat  flux  parameter  for  the  increment  dx  containing  x  =  a.  It 
may  be  more  convenient  for  computation  if  equation  4.29  is  re-written, 
b  . 

q*„ 


w(b) 


•  1 


1 


_d_ 

dx 


P  u  C  vC./2 
1  P  f 


dx 


(4.30) 


In  evaluating  equation  4. 30,  however,  it  must  be  noted  that  if  there  is  a  step  in 

q*  /pu.C  ^Cj2 ,  this  step  contributes  an  increment  of  T  ...  given  by 
nw  ’  p  f  r  w(b)  6  J 

q,: 


AT 


w(b)  ( 


St 


KCJ2 


)  “  (— 
/_  u  Vpu, 


J 


(4.31) 


lCpVc72 

finite  incr 

(b)  Determination  of  temperature  within  the  boundary  layer. 


a,  b  '  ’  1  p  f’  'a 

In  equation  4.31  the  symbol  A  signifies  a  finite  increment. 


The  previous  section  showed  the  method  of  computing  wall  temperature 
consequent  on  an  arbitrary  distribution  of  wall  heat  flux.  A  further  problem  iB, 
however,  the  determination  of  temperature  within  the  boundary  layer.  Again 
the  method  will  be  seen  more  easily  if  the  computation  for  step  q^/pu,C  mZ^/2 
is  first  considered.  ^ 


Temperatures  within  the  boundary  layer  have  been  presented  in  Figs.  4.1, 

4.2,  and  4.3  as  0(x+,  ?,  Pr).  6,  defined  as  in  the  list  of  notation,  is  in  fact 

temperature  on  a  scale  such  that  mainstream  temperature  is  zero  and  wall 

temperature  is  unity.  The  temperature  at  a  point  in  the  boundary  layer  is 

defined  if  T  and  6  be  known.  The  relation  is  T  =  6T  . 

w  w 


Concentrating  our  attention  bn  the  case  of  constant  q^  lpuC^C^/2,  the 

temperature  difference  between  wall  and  mainstream  at  x  =  b,  called  T^^j, 

may  be  found  by  the  method  of  the  previous  section.  The  problem  now  is  to 
find  the  fluid  temperature  at  point  c  with  co-ordinates  x,  y.  Of  the  parameters 
in  Figs.  4.1,  4.2  and  4.3,  x+  and  Pr  are  already  known.  To  determine  6 
and  thus  permit  the  determination  of  temperature  at  y,  the  dimensionless  co¬ 
ordinate  ?  must  be  determined  for  the  point  c.  £  is  defined  by  equation  4.  8  and 
equations  4.8,  4.7,  4.6,  4.5,  4.12,  4.  13  and  4. 14  permit  the  function  £(y+.  Pr) 
to  be  computed.  Results  of  this  computation  are  shown  in  Fig.  4.8  and  Table  4.11. 
Using  Fig,  4.8,  ?  may  therefore  be  found,  after  y+  has  been  determined  from 
equation  4.4.  Thence  6  may  be  read  off  Figs.  4.1,  4.2  or  4.3.  The  difference 
of  temperature  between  point  c  and  the  mainstream  is  then 


-(b) 


ya,  b ,  c) 


a,b 


) 

x  pu.C  ^C„/2 ' 
P  f 


(4.32) 
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In  equation  4.32,  8^  ^  ^  is  the  value  of  6  at  x+  =  x^  ^  and  the  appropriate  ?. 

(St/^  C./2)  ,  is  the  value  read  off  Fig.  4.4  at  x+  =  x+  and  q"/pu,C  V" C  / 2 
f  a,b  a,b  w  «  p  f 

is  the  constant  value  of  the  heat  flux  parameter  between  x  =  a  and  x  =  b. 


When  there  is  an  arbitrary  distribution  q  (x)  the  distribution  of  the  heat 

J  w 

flux  parameter  may  be  computed  as  described  in  the  previous  section  of  the  paper, 
and  the  value  of  the  difference  in  temperature  between  point  c  and  the  mainstream 
computed  by  the  Stieltjes  integral: 


T 


c(b) 


b 

f* 

1 

e 

(a, b,  c) 

d 

.  U 

/  L,w 

I 

o 

(  f  )a,b 

(4.  33) 


Equation  4.  33  is  the  general  equation  by  which  the  results  of  this  paper  may 
be  applied,  since  it  comprehends  equation  4.29  in  that  at  the  wall,  6  is  unity. 
Again  it  may  be  more  convenient  for  computation  if  the  equation  is  rewritten 


'  e(b) 


b 

/ 

°(a,b,e)  d 

•  H 

r  qw  i 

1 

O 

i  ' Lb 

*  u,cpv'r72 

dx 


(4.34) 


In  the  evaluation  of  equation  4.  34  care  must  be  taken  in  allowing  for  steps  in 


<  /*u,c p^V2 


as  was  pointed  out  after  equation  4.30. 


+ 

X 

St  i 

Cx> 

1  t)  u+  a*  d£ 

J  o 

%  diff . 

2 

10 

3.810 

3.838 

+  0.72 

10* 

17.97 

18.07 

+  0.54 

10s 

97.15 

97.21 

+  0.05 

10* 

653.6 

652.4 

-  0.18 

10* 

4952 

4876 

-  1.54 

10* 

39660 

39780 

+  0.30 

Table  4.5  Integrated  Check  of  the  Hesults  for  Pr  =  0.7 


x+ 

St 

x+ 

St 

VcfT2 

^Cf/2 

0 

1 

0.801678 

1.10* 

0.097158 

2 

0.842647 

2 

0.084421 

3 

0.563886 

3 

0.078568 

4 

0.513691 

4 

0.074916 

5 

0.477754 

5 

0.072325 

6 

0.450215 

6 

0.070346 

7 

0.428144 

7 

0. 068760 

8 

0.409884 

8 

0.067446 

9 

0.394414 

9 

0.066330 

10 

0.381063 

1.10* 

0.065363 

20 

0.300255 

2 

0.059528 

30 

0.263463 

3 

0.056654 

40 

0.240197 

4 

0.054776 

50 

0.223636 

5 

0.053401 

60 

0.211018 

6 

0.052325 

70 

0.200964 

7 

0.051448 

80 

0.192695 

8 

0.050711 

90 

0.185731 

9 

0.050077 

1.10* 

0.179756 

1.105 

0,049523 

2 

0.146030 

I 

2 

0.045895 

3 

0.130328 

3 

0.044175 

4 

0.120738 

4 

0  043020 

5 

0.114091 

5 

0.042159 

6 

0.109125 

6 

0.041476 

1 

0.105228 

7 

0.040913 

8 

0.102059 

8 

0.040436 

9 

0.099413 

9t 

0.040023 

io3 

0.097158 

io6 

0.039660 

Table  4.8  The  Spalding  function  St/vrCj,/2  for  Pr  =  0.7 


0.61424 

0.49137 

0.43359 

0.39639 

0.36957 

0.34889 

0.33224 

0.31841 

0.30667 


0.29650 

0.23723 

0.20810 

0.18965 

0.17709 

0.16652 

0.15856 

0.15201 

0.14650 


0.14177 

0.11513 

0.10281 

0.09510 

0.08998 

0.08621 

0.08327 

0.08091 

0.07894 

0.07728 


0,07728 

0.067892 

0.063685 

0.061115 

0.059300 

0.057914 

0.056803 

0.055880 

0.055095 


0.054413 

0.050335 

0.048244 

0.046863 

0.045843 

0.045042 

0.044385 

0.043830 

0.043351 


0.042931 

0.040295 

0.038927 

0.038007 

0.037319 

0.036774 

0.036324 

0.035943 

0.035612 

0.035321 


Table  4.9  The  Spalding  function  St/  for  Pr  =  1 
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300  the  following  relations  hold: 

f+  =  0.3255  exp.  0.4098  ?  (Pr  =  0.7) 

,r+  =  0.09177  exp.  0.4093  ?  (Pr  =  1.0) 

ir+  =  6.317  x  10‘9  exp.  0.4088  ?  (Pr  =  7.0) 
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5.  Description  of  Apparatus 
5.1.  Introduction 


The  purpose  of  this  experiment  was  to  obtain  an  experimental  confirmation  of 
the  numerical  solution  that  is  obtained  for  an  arbitrary  heat  flux  problem.  To 
obtain  this  an  experimental  rig  was  designed  in  which  carbon  dioxide  was  injected 
into  a  pipe  from  a  porous  wall,  the  main  stream  of  the  pipe  being  air.  The  mixed 
gases  were  sampled  for  concentration  of  carbon  dioxide  at  various  stations,  and 
analysed  by  means  of  an  infra-red  gas  analyser;  the  velocity  profile  in  the  pipe, 
just  before  the  injecting  section,  was  fully  developed.  This  was  checked  by 
traversing  a  probe.  (Figs.  5. 2  and  5.  7). 

The  general  layout  of  the  experimental  rig  is  shown  in  Fig.  5.1  and  5.4.  The 
complete  apparatus  can  be  divided  into  the  following  five  groups: - 

1.  Air  section 

2.  Diffusing  section 

3.  Test  sections 

4.  CC>2  section 

5.  Measurements 

5.2.  Air  Section 

The  test  rig  was  on  the  suction  side  of  an  Allis  Chalmers  compressor  with  a 
suction  capacity  of  3120  cu.ft/min.  Air  was  sucked  from  the  atmosphere,  and 
the  air  mass  flow  was  controlled  by: 

1 .  A  control  valve  on  the  by-pass  line 

2.  A  control  valve  in  the  rig  after  the  orifice 

Both  the  control  valves  were  butterfly  valves  with  angular  indicators  and  stop  nuts 
on  handle  for  fixing  the  valve  positions. 

The  pipe  size  was  made  4.75"  internal  diameter  for  the  following  two  reasons:  - 

1.  To  facilitate  the  traverse  of  the  probe  for  measuring  concentration 
profile  very  accurately,  it  was  better  to  have  a  fairly  large  size 
pipe  diameter. 

2.  The  final  size  4.75"  I.D.  was  decided  by  the  size  of  the  porous  tube 
available. 

To  obtain  a  fully  developed  and  an  axisymmetric  flow,  a  bell -mouth  entry  was 
chosen.  The  bell-mouth  section  was  made  from  wood  and  polished  to  obtain  a  very 
smooth  surface. 

There  is  a  possibility  of  asymmetric  flow  in  the  pipe  if  there  is  the  slightest 
error  in  coincidence  of  the  axis  of  the  bell-mouth  section  with  that  of  the  main  pipe. 
To  prevent  this,  a  straightening  grid  was  inserted  immediately  after  the  bell -mouth 
section.  The  straightening  grid  was  made  by  inserting  172  stainless  steel  tubes, 
0.323"  O.D.,  0.007"  thick  and  2"  long,  in  a  flange  2"  thick,  4.75"  I.D.  (Fig.  5.6). 
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The  gap  was  filled  by  a  tube  separately  made  to  suit  the  gap.  This  tube  was  0.412" 
O.D.  0.008"  thick  and  2"  long. 

The  inlet  section  was  made  22' -6"  long  and  in  three  pieces,  9  ft.,  6ft.  and 
7ft.  -6"  respectively.  This  long  length  of  inlet  section,  approximately  57  times 
the  internal  diameter  of  the  tube,  was  used  to  ensure  that  the  air  flow  just  before 
the  test  section  and  the  diffusing  section,  was  fulLy  developed  turbulent.  Further, 
to  obtain  a  very  symmetrical  flow,  the  last  piece  of  inlet  section,  7ft. -6"  long,  and 
before  the  diffusing  section,  was  made  from  perfectly  circular  drawn  tube  of  one 
length  and  0.125"  thick.  The  other  two  pieces  were  made  from  16  SWG.  3ft.  wide 
steel  sheet. 

Between  the  test  section  and  orifice,  a  straight  pipe  7ft.  long  was  connected  in 
accordance  with  the  British  Standards  Specification  for  orifices.  (Fig.  5.8). 

5.3.  Diffusing  Section 

The  diffusing  section  mainly  consisted  of  a  porous  tube  and  a  diffusing  chamber. 
The  porous  tube  was  a  high  grade  ceramic  tube  with  uniform  porosity. 

The  aim  of  the  experiment  was  to  measure  the  concentration  profiles  of  the 
diffusing  substance  at  different  sections  for  uniform  step  mass  transfer,  followed 
by  a  step  adiabatic  wall.  Mass  transfer  through  a  porous  material  depends  on 
three  factors: - 

1 .  Static  pressure  difference 

2.  Thickness  of  the  porous  tube 

3.  Porosity  of  the  material 

Hence  to  obtain  uniform, 

1.  pressure  in  the  chamber, 

2.  thickness  of  the  material,  and 

3.  porosity  throughout  the  diffusing  length, 

the  following  construction  was  used:- 

1.  High  grade  ceramic  tube  was  used  to  ensure  uniform  porosity. 

2.  The  outside  and  inside  surfaces  were  ground  to  uniform  diameters. 
During  grinding  precautions  were  taken  to  obtain  both  the  surfaces 
co-axial  and  hence  of  uniform  thickness.  The  manufacturers  of  the 
ceramic  tube  confirmed  that  grinding  does  not  alter  the  porosity. 

3.  The  diffusing  chamber  was  made  12"  l.D.  and  10.25"  wide.  This 
large  size  of  chamber  gave  enough  room  for  the  carbon  dioxide  to 
settle  before  it  diffused  through  the  porous  tube.  This,  in  other 
words,  ensured  uniform  pressure  in  the  chamber. 

4.  Carbon  dioxide  was  admitted  into  the  chamber  from  16  places,  eight 
on  either  side  of  the  chamber,  equiplaced  on  11.125"  P.C.D.  Two 
circular  ring-mains,  made  from  0.5"  diameter  copper  tube,  supplied 
carbon  dioxide  to  the  chamber  at  these  16  positions  by  0.25"  dia. 
tubes  as  shown  in  Fig.  5.5.  Carbon  dioxide  to  the  ring-mainB  was 
supplied  from  a  carbon  dioxide  cylinder,  through  the  pressure 
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5.4. 


reducing  valve,  main  valve  and  the  flowmeters. 

5.  'O'  ring  seals  were  used  at  all  points  to  prevent  any  leakages. 

Test  Sections 


There  were  two  test  sections,  a  small  one  for  measuring  the  velocity  profile 
and  one  large  one  for  measuring  the  concentration  profiles.  They  were  made  from 
4.5"  I.D. ,  0.281"  thick  mild  steel  tubes,  and  they  were  bored  accurately  to 
4.75"  internal  diameter.  The  flanges  of  these  test  sections  were  made  to  fit  in  a 
circular  groove  of  the  adjoining  flanges  so  that  they  could  be  rotated  around  the 
axis  of  the  pipe.  With  this  arrangement,  traversing  could  be  made  at  different 
circumferential  positions  to  check  the  axisymmetry  of  the  flow. 

The  smaller  section  was  mounted  before  the  diffusing  chamber.  It  was  made 
8"  long  with  only  one  boss  welded  in  the  middle  of  the  length  to  take  the  traversing 
assembly. 

The  larger  section  was  mounted  after  the  diffusing  chamber.  It  was  made 
19"  long  with  four  bosses  welded  at  four  axial  positions  and  at  right  angles  in  a 
consecutive  order  (see  Fig.  5.5). 

During  the  construction  of  these  test  sections,  all  care  and  precautions  were 
taken  to  ensure  that 

1.  The  axes  of  these  test  sections  were  in  line  with  the  axis  of  the 
main  pipe  line. 

2.  No  leakages  occurred  from  side  surfaces  where  they  rotate.  To 
ensure  this  'O'  ring  seals  were  used  at  thesesurfaces. 

3.  Traversing  axis  of  the  probe  was  perpendicular  to  the  axis  of  the 
pipe  and  in  the  same  plane . 

4.  Surfaces  of  the  matching  pieces,  on  which  the  traversing  probe 
assembly  was  mounted,  were  parallel  and  at  a  fixed  distance 
from  the  axis  of  the  pipe. 

5.  Inside  surfaces  of  the  matching  pieces  on  which  the  traversing  probe 
assembly  was  mounted,  matched  the  inside  curvature  of  the  pipe. 

To  obtain  this,  test  sections  were  bored  accurately  to  4.75"  I.D.  , 
with  matching  pieces  in  their  position. 

6.  Inside  surfaces  was  smooth  and  straight. 

5.5.  CC>2  Section 

CO2  was  obtained  from  three  cylinders  (each  containing  28  lb  of  liquid  CO2) 
connected  in  parallel.  From  these  cylinders,  CO2  passed  through: - 

1.  Pressure  regulating  valve 

2.  Settling  chamber 

3.  Control  valve 

4.  Two  flowmeters ,  and 

5.  Orifice 

before  reaching  the  diffusing  section,  ^''copper  tubes  were  used  for  CO2  pipe  line. 
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At  the  pressure  reducing  valve,  the  temperature  of  the  CO2  was  very  low 
due  to  its  expansion.  To  bring  the  temperature  of  CO2  back  to  atmospheric 
temperature  at  flowmeter  and  at  diffusing  sections,  the  CO2  cylinders  and 
pressure  reducing  valve  were  placed  at  a  distance  from  the  test  rig  and  the  CO2 
cylinders  were  placed  in  a  water  tank.  With  this  arrangement,  CO2  passed 
through  the  long  length  of  copper  tube  (exposed  to  atmosphere)  and  in  consequence 
obtained  atmospheric  temperature  before  reaching  the  test  rig. 

5.6.  Measurements 


(a)  Air  mass  flow;  Air  mass  flow  was  measured  by  an  orifice  at  the 
downstream  side  of  the  test  section.  The  orifice  plate  was  of  "D  and  D/2  taps" 
type,  and  it  was  made  according  to  the  British  Standard  Specification  B.S.  1042 
(1943).  To  simplify  the  calculations,  the  area  ratio  of  orifice  to  pipe  was  made 
0.5.  There  were  four  sets  of  pressure  tappings  at  four  right  angle  circumferential 
positions.  Pressure  differences  across  the  orifice  were  measured  on  vertical 
water  manometers.  The  mean  of  these  four  sets  of  pressure  differences  was 
used  to  obtain  air  mass  flow. 

(b)  CO2  flow:  The  flow  of  CC>2  was  measured  by  two  Fischer  &  Porter 

percentage  type  flowmeters.  They  read  accurately  within  t  2%  of  the  maximum 
flow.  Two  flowmeters  were  used  in  series  to  cover  the  wide  range  of  CO2  flow. 
100%  reading  on  the  small  flowmeter  represented  0.97  cu.ft.  /min.  at  14.7  p.s.  i. 
and  60°F.  ,  and  100%  reading  on  large  flowmeter  represented  3.72  cu.ft. /min. 

at  14.7  p.s.i.  and  70°F.  For  different  temperature  and  pressure  conditions,  the 
flow  was  corrected  by  the  use  of  curves  supplied  by  the  manufacturers. 

At  first  it  was  found  that  the  amount  of  CO2  flowing  through  the  pipe  and 
measured  by  the  gas  analyser  was  less  than  the  CO2  measured  by  the  flowmeter. 
After  checking  leaks  in  the  CO2  line,  the  flowmeters  were  suspected.  For  an 
approximate  check  on  flowmeters,  an  orifice  was  introduced  in  the  CO2  pipe  line, 
and  this  confirmed  the  accuracy  of  the  flowmeters  (see  Figs.  7.3  and  7.4).  Later 
on  the  fault  was  found  in  the  sampling  cell  of  the  infra-red  gas  analyser.  The  cell 
was  cracked  and,  in  consequence,  air  was  leaking  into  the  cell  resulting  in  low 
concentration  measurements.  The  cell  was  then  replaced  by  a  new  one. 

(c)  Concentration  of  CO2:  Concentration  of  CO2  was  measured  by  an 

infra-red  gas  anAlySer'Ty'pe  S.B.  1 .  This  gas  analyser  has  two  absorbing  lengths. 
One  length  is  a  standard  cell,  whereas  the  other  length  has  two  cells,  one  for  pure 
air  and  another  for  the  sample  to.be  tested.  Radiation  from  the  nichrome  heaters 
passes  through  the  absorption  tubes  and  thence  two  two  receivers,  which  are  filled 
with  the  gas  to  be  detected.  These  receivers  are  partitioned  off  from  one  another 
by  a  thin  metal  diaphragm,  and  form  an  electric  condenser.  Absorption  of  infra¬ 
red  radiation  in  the  sampling  cell  depends  on  the  concentration  of  carbon  dioxide. 

This  absorption  of  radiation  causes  a  pressure  difference  between  the  two 
receiving  chambers,  with  the  result  that  the  thin  diaphragm  deforms  and  causes 
a  variation  of  the  capacity.  The  resulting  capacity  changes  are  amplified 
electronically  and  a  final  indication  is  obtained  on  an  output  galvanometer. 


The  sampling  length  was  on  the  suction  side  of  a  Charles  Austen  Mk.  1  pump 
with  a  capacity  of  3500  c.c./min.  The  amount  of  suction  was  adjusted  by  controlling 
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the  speed  and  the  opening  on  the  by-pass  line.  The  sample  from  the  traversing 
probe  was  passed  through  the  dr  ying  tube  and  filter  before  passing  through  the 
sampling  cell.  Pure  air  was  passed  through  the  drying  tube  and  a  tube  containing 
soda  asbestos  to  remove  carbon  dioxide  before  it  entered  the  standard  cell 
(Fig.  5.3).  Moisture  was  removed  from  the  sample  because  it  absorbs  the  infra¬ 
red  radiation  of  the  same  wavelength  as  carbon  dioxide  and  would  give  a  false 
reading.  Thin  (1/8"  O.D.)  polythene  tubing  was  used  for  connections. 

The  gas  analyser  was  calibrated  for  two  ranges  of  concentration,  and  the 
calibration  curves  are  shown  in  Figs,  7.  1  and  7.2.  Calibration  was  done  by 
mixing  a  known  amount  of  carbon  dioxide  with  a  known  amount  of  air,  and  measuring 
the  galvanometer  reading  for  the  mixture.  The  flow  of  carbon  dioxide  was  measured 
by  a  soap-bubble  flowmeter  and  airflow  was  measured  by  a  gasmeter. 

During  calibration  a  standard  bottle  was  made  in  which  carbon  dioxide  was 
mixed  with  air,  and  the  galvanometer  reading  for  that  mixture  was  measured. 

The  instrument  was  also  provided  with  a  calibration  system  consisting  of  a  metallic 
wire  which  comes  between  the  sampling  tell  and  a  receiving  chamber  (when 
operated  by  a  solenoid  switch!  and  absorbs  a  fixed  amount  of  radiation  giving  a 
standard  galvanometer  reading. 

During  tests,  the  same  conditions  and  suction  were  used  as  during  calibration, 
the  instrument  being  frequently  checked  by  the  "wire"  reading  and  the  standard 
bottle  for  the  same  galvanometer  readings. 

(d)  Pressure  measurement:  To  obtain  the  friction  factor  in  the  pipe,  the 

static  pressure  gradient  for  a  fixed  length  of  pipe  was  measured  on  a  Betz  mano¬ 
meter.  This  length  was  selected  where  the  flow  was  fully  developed  (after  a  length 
of  40  diameters  from  entry).  The  velocity  profile  was  measured  by  measuring 
the  static  pressure  on  a  vertical  water  manometer  and  the  difference  of  pressure 
between  the  total  and  static  on  a  Betz  manometer.  The  pitot  tube  was  traversed 
at  the  test  section  before  the  diffusing  chamber. 

6.  Test  Procedure  -  Calculation  Procedure 

8.1.  Introduction 

The  aim  of  this  experimental  investigation  was  to  obtain  an  experimental 
confirmation  of  the  numerical  solution.  The  test  rig  was  made  in  which  carbon 
dioxide  was  injected  through  a  porous  section  of  the  tube  wall,  while  air  was 
flowing  turbulently  in  the  axial  direction.  The  injecting  length  was  10.25"  and 
the  concentration  profiles  in  the  boundary  layer  were  measured  downstream  of 
the  injection  section.  The  injection  was  assumed  to  be  uniform.  Conditions  are 
therefore  such  that  the  concentration  profiles  are  measured  for  a  step  length  of 
uniform  injection  followed  by  an  impermeable  wall.  In  heat  transfer,  this  is  a 
problem  of  finding  temperature  profiles  for  step  heat  flux  followed  by  an  adiabatic 
wall.  The  flow  in  the  pipe  was  fully  developed  and  wall  shear  stress,  therefore,  was 
uniform  in  the  test  section. 

8.2.  Procedure  of  Test 

The  test  procedure  was  very  simple.  It  involved  traversing  a  probe  and 
measuring  the  galvanometer  reading,  which  in  consequence  gave,  directly  from 
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the  calibration  curves,  the  concentration  of  the  carbon -dioxide  at  that  particular 
position,  A  test  card  (Fig.  6.0)  is  attached  to  explain  the  measurements  taken 
during  the  tost.  Each  test  took  approximately  two  days  to  run,  and  comprised: 

1 .  Three  to  and  fro  traverses  for  concentration  profiles  at  the 
axial  section  but  at  three  circumferential  positions  called 
positions  A,  B,  and  C. 

2.  Velocity  traverse  before  the  diffusing  section. 

During  the  test,  the  main  stream  air  flow  and  the  flow  of  carbon  dioxide 
were  kept  constant.  All  possible  precautions  were  taken  to  obtain  true  readings, 
e.g. 


1 .  Giving  enough  time  for  the  galvanometer  reading  to  Bettle 

2.  Checking  drying  agents 

3.  Checking  very  often  the  wire  reading  of  the  I.R.G.A. 

4.  Checking  air  flow  and  carbon  dioxide  flow,  etc. 

Before  starting  the  actual  profile  tests,  many  preliminary  tests  were  carried 
out  to  ensure  that  all  the  measurements  were  correct  and  reliable. 

6.3.  Preliminary  tests 

To  confirm  that  the  flow  was  fully  turbulent,  the  velocity  profiles  for  various 
main  stream  velocities  were  measured  at  the  test  section  before  the  diffusing 
section.  The  results  were  found  to  lie  on  a  plot  of  Universal  Velocity  Profile 
(Fig.  7.5).  Velocity  profiles  obtained  were  integrated  to  check  the  air  mass 
flows,  and  were  found  to  agree  within  2%  with  the  mass  flows  obtained  from  the 
orifice  meter.  Simultaneously,  ratios  of  average  to  maximum  velocities  were 
calculated  and  compared  with  the  previous  investigators.  The  results  were  found 
to  be  satisfactory  (Fig.  7.8). 

The  knowledge  of  the  friction  factor  was  necessary  for  the  theoretical 
calculation  of  the  concentration  profile.  Instead  of  using  the  well  known  equation 
of  Von  Karman  etc.  for  the  calculation  of  the  friction  factor,  it  was  decided  to 
measure  the  actual  friction  factor  and  use  this  experimental  value  for  the 
theoretical  calculation.  The  friction  factor  was  obtained  by  measuring  the 
pressure  drop  in  the  pipe, across  a  length  of  57",  on  the  Betz  manometer.  The 
results  are  shown  in  Fig.  7.9.  The  best  fit  curve  through  the  experimental 
points  was  found  to  be 

Cf/2  =  0.04135  (Red)'0'254  (6.1) 

The  gas  analyser  was  calibrated  for  two  ranges  of  concentration  0.  to  1.1% 
and  0  to  2.9%.  To  prevent  any  error,  the  calibrations  were  carried  out  in  the 
conditions  that  would  occur  during  the  concentration  profile  test.  Calibration 
was  carried  out  by  the  method  indicated  earlier. 

To  measure  the  flow  of  carbon  dioxide,  two  flowmeters  were  used  in  series 
to  cover  a  wide  range.  These  flowmeters  were  compared  with  each  other  in  the 
intermediate  region  of  flows  for  which  both  gave  indications  (Fig.  7.3). 
Simultaneously,  an  orifice  was  included  in  the  line  to  check  the  flow  measured  by 
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the  flowmeter.  The  orifice  was  not  made  according  to  the  British  Standard 
Specification  which  limits  to  the  orifice  diameter  of  0.5".  In  this  instance  the 
diameter  of  the  orifice  was  0.25"  with  m  =  0.346.  For  such  a  small  diameter 
the  correction  factor  was  very  high  and  was  out  of  the  region  given  in  B.S.  1042. 
The  flowmeter  therefore  was  compared  with  the  orifice  meter  without  applying  the 
correction  factor  to  the  orifice  readings.  The  error  in  the  orifice  meter  was 
therefore  expected  to  be  in  the  region  of  10  to  15%.  Comparison  is  shown  in 
Fig.  7.4. 


It  is  known  that  the  concentration  at  any  point  in  the  boundary  layer  is 
directly  proportional  to  the  injection  rate.  To  check  that  the  measurements  are 
reliable,  a  quick  test  was  carried  out  in  which  the  concentrations  were  measured 
at  a  fixed  position  in  the  boundary  layer  for  various  injection  rates.  The  test 
was  carried  out  for  two  different  main  stream  velocities.  The  results  confirmed 
the  reliability  of  the  measurements  (Fig.  7.7.). 

6.4.  Calculation  Procedure  (Experimental) 

The  concentration  of  the  injected  carbon  dioxide  for  a  given  galvanometer 
reading  was  obtained  directly  from  the  calibration  curves.  For  each  test  there 
were  three  sets  of  measurements  for  concentration  profile  at  three  circumferential 
positions,  and  a  mean  curve  through  these  three  sets  of  results  was  considered 
as  a  concentration  profile  at  that  particular  section  for  that  particular  main  stream 
velocity  and  injection  rate. 


The  average  velocity  of  air  through  the  pipe  was  calculated  from 

(6.2) 

h  =  hs  -  h6  =  pressure  differential  across  the  orifice  in  inches  of  water 
Ts  =  temperature  of  air  at  upstream  of  the  orifice  in  °K. 

P  =  pressure  upstream  of  the  orifice  in  lb./in^. 

Z  =  correction  factor  from  B.S.  1042. 


Reynolds  number  was  obtained  from 


Red 


uD 

v 


u  x  4.75 
12  v 


(6.3) 


Flow  of  carbon  dioxide  V  in  cu.ft.  /sec.  was  obtained  by  multiplying  the 
percentage  reading  on  the  flowmeter  with  a  constant  factor.  100  per  cent  reading 
on  the  large  flowmeter  represented  0.062  cu.ft. /sec.  at  14.7  p.s.i.  and  70  F. 
and  100  percent  reading  on  the  small  flowmeter  represented  0.0162  cu.ft. /sec. 
at  14  7p,s.i.  and  70  F.  Correction  factors  for  temperature  and  pressure  were 
obtained  from  the  curves  supplied  by  the  manufacturers.  The  injecting  surface 
area  of  the  porous  tube  was  1.062  ft  . 

V  "ifc-T  cu.  ft. /sec.  ft2.  (6.4) 

w  1.062 


To  compare  the  experimental  concentration  profiles  with  the  theoretical 
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concentration  profiles,  the  readings  of  the  concentration  profiles  were  multiplied 

by  a  factor  — —  .  M  .  In  other  words,  all  the  concentration  profiles  were  for 
100  Vw 

a  ratio  of  injection  velocity  to  main  stream  velocity  =  1/100. 

The  velocity  profile  in  the  boundary  layer  was  obtained  by  traversing  a  probe 

and  using  _ 

I  hT 

u  =  14.93  J  ~  (6.5) 

where  h  =  h3  -  h^=  total  pressure  -  static  pressure  in  inches  of  water 

T„  =  temperature  in  °K 
,  2 

P2  =  Pressure  in  lb. /in  . 

During  the  friction  factor  test,  the  friction  factor  was  calculated  from 


l.  e. 


—  x  t/DL  =  (P  -  PJ  ~  D2 
g  124 

C/2  -  81  P<  ~  P~  ’  *  D. 

L  =  length  of  pipe  in  ft. 


D  =  diameter  of  the  pipe  in  ft. 

2 

(F|  -  P2)  =  pressure  loss  in  lb. /ft  . 


(6.6) 

(6.7) 


6.5. 


Calculation  Procedure  (Theoretical) 


Section  1 . 


Where  concentration 


Fig.  6.1 
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Fig.  6. 3 


In  the  experimental  rig,  the  injection  was  assumed  to  he  uniform  and 
therefore  this  was  a  case  of  step  heat  flux  followed  by  an  adiabatic  wall 
(Fig.  6.1).  To  apply  the  numerical  solution  for  step  heat  flux,  this  problem 
was  considered  in  two  steps: 

1.  Solution  for  uniform  injection  from  A  to  C  (Fig.  6. 2) 

2.  Solution  for  uniform  injection  from  B  to  C  (Fig.  6.3) 

The  final  solution  was  obtained  by  subtracting  the  second  solution  from  the 
first.  For  carbon  dioxide  diffusing  in  air,  the  Schmidt  number  is  0.91,  and 
therefore  the  numerical  solution  for  Pr  =  1  was  used  to  obtain  the  theoretical 
concentration  profile. 

The  calculation  of  the  theoretical  concentration  profile  can  be  summarised 
as  follows:  - 

1 .  Calculation  of  x|  and  x  *  from 

x+  =  x  Red  x  5  (6.8) 

L  =  injecting  length. 

u  and  Red  are  known  from  the  experimental  calculation.  Friction  factor  C^/2 
is  obtained  from  the  experimental  relation: 

Cf/2  =  0.04135  (Red)' 


(6.9) 
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(2l  Reading  the  values  of  Spalding  functions: 


St 


,  and 


for  x  and  x2  from  Fig.  4.4. 


(31  Knowing  ^C/2 ,  (St  )  and  (St  )  are  obtained  from  (2). 
f  ml  m2 

V  " 

W  J|« 

St  =  -  =  Equivalent  of  Stanton  number  in  mass  transfer 

m  - 

C  u 
w 

,  2 
V  "  =  iniection  rate  in  cu.  ft. /see.  ft  . 
w  J 

u  =  average  velocity  ft. /see. 

C  =  surface  concentration  of  the  carbon  dioxide  in  cu.  ft. /cu.  it. 
w  ,  •  . 

of  mixture. 

V  " 

(4)  Theoretical  solution  was  carried  out  for  =  1/100.  Hence 

u 

C  and  C  are  calculated  from  (3). 

W1  W2 

(5)  For  Schmidt  number  equal  to  1,  the  concentration  profiles  are  given 
in  Fig.  4,  in  the  form  of  C/C'w  for  a  large  range  of  x4  .  Knowing 

and  C  ,  the  concentration  profiles  for  x4  and  x4  are  obtained. 

W2  r  12 

(6)  The  final  theoretical  concentration  profile  is  obtained  by  subtracting  the 
concentration  profile  for  x4  from  the  concentration  profile  of  xf  . 

(7)  The  concentration  profile  obtained  above  is  in  the  form  of  C(?)  and 

not  C(y).  From  Fig.  4.8  and 

y+  =  o  ~  <6.10) 

C(y)  can  be  obtained. 


Note  that  this  definition  is  identical  with  St  -  - =r  where  m  is  the 

m  p  u  w 

w 

mass  rate  of  flow  of  carbon  dioxide  per  unit  area  at  the  wall,  and  p  is  the 

v  w 

true  density  of  carbon  dioxide  at  the  wall  (i.e.  mass/unit  volume  of  space). 
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7.  Results  and  Conclusions 
7.1.  Results 


The  results  of  the  preliminary  tests  and  eight  concentration  profile  tests 
are  shown  in  fourteen  tables  (Tables  7.1  to  7. 14)  and  21  graphs  (Figs.  7. 1 
to  7. 21). 

Tables  7.1  to  7.8  contain  the  measurements  of  eight  concentration  profile 
tests.  Table  7.  9  contains  the  summary  of  the  eight  tests.  Measurements  for 
the  calibration  of  the  I.  R .  G.  A .  are  shown  in  Tables  7.10  and  7.11.  Tab  le  7.12 
shows  the  measurements  of  friction  factor,  and  Tables  7.13  and  7.14  show  the 
comparison  of  the  flowmeters  for  CC^. 

The  calibration  of  the  I.  H.  G.  A.  is  shown  in  Figs.  7.1  and  7.2.  Figs.  7.3  and 
7.4  show  comparison  of  flowmeters  for  carbon  dioxide.  Measurement  of  the 
velocity  profile  is  shown  in  Fig.  7.5,  and  comparison  with  Spalding's  velocity 
profile  is  given  in  Fig.  7.6.  The  results  of  the  friction  factor  measurements  are 
shown  in  Fig.  7.9.  In  Fig.  7.8,  the  ratio  of  the  average  to  maximum  velocity 
is  compared  with  the  results  of  the  previous  investigators.  The  results  of  the 
test  confirming  the  direct  variation  of  the  concentration  with  the  injection  rate 
is  given  in  Fig.  7.7.  In  Fig.  7.10,  the  theoretical  concentration  profile  at 
section  1  is  shown  for  various  main  stream  velocities. 

Figs.  7.11  to  7.  18  show  the  experimental  concentration  profiles  of  the  eight 
tests.  Fig.  7J.9  show  the  concentration  profiles  for  tests  5  to  8,  in  which  the 
main  stream  velocity  was  fixed  but  the  injection  rate  was  varied,*  Fig.  7.20 
gives  the  concentration  profiles  of  tests  5  to  8  (converted  for  fixed  V'J’  )  confirming 

the  linearity  of  the  concentration  with  the  injection  rate.  Finally,  Fig.  7.21  shows 
the  comparison  of  the  experimental  concentration  profiles  with  the  theoretical  profiles. 

7.2.  Conclusions 


On  the  theoretical  side,  the  problem  of  heat  and  mass  transfer  from  non- 
isothermal  surfaces  has  been  attacked  in  two  ways. 


1.  Approximate  solutions  are  obtained  for  laminar  and  turbulent  cases,  by 
assuming  the  shape  of  velocity  and  temperature  profiles  in  the  boundary  layer 
and  integrating  the  energy  equation.  The  equations  obtained  are  very  simple  and 
predict  results  that  are  in  agreement  with  those  of  Rubesin'^  ,  and  Reynolds  et 
al(10V 


2.  Exact  solution  of  the  partial  differential  equation  is  obtained  for  a  step  heat 
flux  problem  and  for  Prandtl  numbers  0,7,  1  and  7.  The  numerical  integration 
was  carried  out  on  a  Digital  Computer.  Spald'ng's  Velocity  Law  was  used  with 
an  assumption  of  unity  turbulent  Prandtl  number.  A  method  has  been  developed 
to  apply  this  solution  to  an  arbitrary  heat  flux  problem.  Various  methods  are 
known  to  exist  in  which  the  heat  transfer  coefficient  can  be  calculated  for  an 
arbitrary  heat  flux  case.  But  so  far  as  the  authors  know,  the  method  presented 
in  this  report  is  the  only  method  that  exists  in  which,  in  addition  to  the  calculations 
of  heat  transfer  at  the  surface,  the  temperature  within  the  boundary  layer  can  be 
calculated.  Due  to  lack  of  time,  it  has  not  been  possible  to  compare  these 
solutions  with  the  experimental  results  of  Reynolds  et  al'10'. 
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TABLE  7.3.  RESULTS  OF  TEST  3  TABLE  7.4.  RESULTS  OF  TEST  4 
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TABLE  7.5.  RESULTS  OF  TEST  5  TABLE  7.6.  RESULTS  OF  TEST  6 
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TABLE  7. 7.  RESULTS  OF  TEST  7  TABLE  7.8.  RESULTS  OF  TEST  8 


Small  Flowmeter 
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TABLE  7.10.  CALIBRATION  OF  I.R.G.A.  (Range  0  -  1.1%) 
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TABLE  7.11.  CALIBRATION  OF  I.R.G.A.  (Range 
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No. 

u 

ft/sec. 

P 

Ibs/cu.  ft. 

Red 

x  10-5 

hf 

in  m .  m . 

Cf/2 

1 

87.7 

0.0753 

2.205 

7.735 

0.001827 

D 

82.4 

0.0754 

2.073 

6.85 

0.001827 

H 

78.1 

0.0755 

1.963 

6.26 

0.001855 

B 

92.1 

0.0752 

2.314 

8.4 

0.0018 

5 

118.9 

0.07435 

2.99 

13.1 

0.001707 

6 

114.7 

0.07445 

2.882 

12.135 

0.001692 

7 

110.4 

0.0746 

2.777 

11.325 

0.001705 

8 

146.5 

0.07335 

3.685 

18.4 

0.001596 

9 

139.7 

0.0736 

3.51 

16.85 

0.001606 

10 

134.7 

0.0738 

3.387 

15.875 

0.001618 

11 

127.9 

0.074 

3.22 

14.5 

0.00164 

12 

102.8 

0.0748 

2.587 

9.85 

0.001705 

13 

73.35 

0.0756 

1.844 

5.2 

0.00175 

14 

60.3 

0.0758 

1.517 

3.975 

0.001974 

15 

49.4 

0.076 

1.243 

2.84 

0.002092 

16 

34.9 

0.0762 

0.8775 

1.58 

0.00232 

17 

24 

0.07625 

0.6035 

0.86 

0.002675 

TABLE  7.12,  FRICTION  FACTOR  MEASUREMENTS 

„  g<P«  '  P.)  D 

1/ 2  =  -a  4L 

P  u 


L 

D 


57" 

4.75" 
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No. 

Reading  on  small 
flowmeter 

Reading  on  large 
flowmeter 

Flow  in  cu.  ft/ mt. 
small  flowmeter 

Flow  in  cu.ft/nr 
large  flowmeter 

1 

93 

24.6 

0.911 

0.915 

I'S 

96 

25.4 

0.941 

0.944 

It 

90 

23.8 

0.882 

0.885 

n 

83.7 

22.2 

0.82 

0.826 

5 

81 

21.6 

0.7935 

0.8035 

6 

75 

19.8 

0.735 

0.7365 

7 

59.5 

15.3 

0.583 

0.569 

8 

97 

25.5 

0.95 

0.948 

9 

89 

23.4 

0.872 

0.87 

10 

87 

23 

0.852 

0.856 

11 

73 

19.2 

0.715 

0.714 

12 

82 

21.8 

0.803 

0.811 

13 

C2.5 

22 

0.808 

0.8185 

14 

79 

20.85 

0.7735 

0.775 

15 

74.6 

19.5 

0.730 

0.725 

16 

68 

17.9 

0.660 

0.666 

■ 

85 

22.4 

0.833 

0.833 

TABLE  7.13.  COMPARISON  OF  FLOWMETERS 


100  percent  reading  on  small  flowmeter  represents 
0.98  cu.  ft/mt  at  14.7  p.  s.  i.  and  70°F. 

100  percent  reading  on  large  flowmeter  represents 
3.72  cu.  ft/mt  at  14.7  p.s.i.  and  70°F. 
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The  experiment  was  started  with  the  idea  of  obtaining  an  experimental 
confirmation  of  the  numerical  solution.  It  was,  therefore,  planned  to  measure 
the  concentration  profiles  at  :  - 

1.  Various  sections  downstream  of  the  injection  for  various, 

2.  mainstream  velocities, 

3.  injection  rates, 

4.  injecting  lengths. 

Eight  tests  were  carried  out  in  which  concentration  profiles  were  measured  at 
section  1  for  a  few  main  stream  velocities  and  injection  rates. 

Comparison  of  the  experimental  results  with  the  theoretical  solution  is 
shown  in  Fig.  7.20  and  Fig.  7.21.  It  is  seen  from  this  comparison  that 

(1)  Experimental  results  show  higher  diffusion  of  CO2  in  the  boundary  layer 
than  that  of  the  theoretical  solution,  and  in  consequence: 

(2)  the  measured  surface  concentration  is  lower  than  that  of  the  theoretical 
solution. 

(3)  The  shape  of  the  measured  concentration  profile  is  different  from  that  of 
the  theoretical  solution. 

There  is  no  doubt  about  the  accuracy  of  the  experimental  results  because: - 

1 .  The  amount  of  carbon  dioxide  obtained  by  integrating  the  concentration 
profile  agreed  with  the  amount  of  carbon  dioxide  injected  through  the 
porous  tube. 

2.  The  concentration  at  any  point  varied  directly  with  the  injection  rate. 

3.  The  shape  of  the  concentration  profile  did  not  vary  with  the  injection  rate. 

4.  All  the  eight  tests  consistently  showed  a  similar  shape  of  the  concentration 
profile  in  the  boundary  layer. 

The  authors  feel  that  the  disagreement  with  the  theoretical  solution  are  due 
to  the  following  reasons: - 

(1)  The  theoretical  concentration  profiles  for  comparison  with  the  experimental 
results  were  obtained  on  the  assumption  that  the  injection  through  the  porous 
tube  is  uniform  throughout  the  injecting  length.  There  is  no  experimental 
confirmation  to  this  assumption. 

(2)  The  theoretical  solution  is  for  two  dimensional  flow,  and  measurement  of 
higher  diffusion  may  be  due  to  reduction  in  cross-sectional  area  (due  to 
decrease  in  radius)  with  increase  in  y  (distance  from  surface). 
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diameter  of  the  pipe  ,  ft. 

distance  parallel  to  the  wall  and  from  the  beginning  of  the 
boundary  layer,  ft. 

distance  where  the  heating  starts,  ft. 

distance  normal  to  the  wall,  ft. 

variable  of  integration 

stream  function 

boundary  layer  thickness,  ft. 

momentum  thickness,  ft. 

shear  thickness,  ft. 

thermal  boundary  layer  thickness,  ft. 

enthalpy  flux  thickness,  ft. 

heat  flux  thickness,  ft. 

velocity  in  the  boundary  layer,  ft.  / sec. 

mean  velocity  parallel  to  wall  in  the  pipe,  ft.  / sec. 

main  stream  velocity  parallel  to  wall,  ft.  / sec. 

velocity  normal  to  the  wall,  ft.  / sec. 

temperature  in  the  boundary  layer,  °F. 

o 

main  stream  temperature,  F. 
temperature  of  the  surface,  °F. 

concentration  of  the  injected  gas  in  the  boundary  layer  in 
cu.  ft.  / cu.  ft. 

concentration  of  the  injected  gas  outside  concentration  boundary 
layer  (this  is  zero)  in  cu.  ft.  / cu.  ft. 

concentration  of  the  injected  gas  at  the  surface  in  cu.  ft. /cu.  ft. 
viscosity,  lb. /ft.  sec. 

2 , 

kinematic  viscosity,  ft  /sec. 
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€ 

2 , 

eddy  viscosity,  ft  / sec. 

h 

n 

heat  transfer  coefficient  BTU/sec.ft  .  F. 

k 

thermal  conductivity  BTU/sec.ft.  °F. 

h 

c 

2 

mass  transfer  coefficient  cu.  ft. /sec.  ft. 

p 

3 

density,  lb. /ft  . 

C 

P 

specific  heat  BTU/lb  °F 

D* 

2 

diffusion  coefficient,  ft  /sec. 

€h 

2 

eddy  thermal  diffusivity  ft  /sec. 

€ 

C 

2 

eddy  mass  diffusivity  ft  / sec. 

a 

2 

thermal  diffusivity,  K/p  C  ,  ft  / sec. 

T 

W 

3  2 

shear  stress  at  the  surface,  lb. /ft  /sec  . 

•  * 

heat  transfer  rate  at  the  wall  per  unit  area  and  unit  time 

•  M 

m 

w 

2 

mass  transfer  rate  at  the  wall  in  lb.  ft/ sec.  ft  . 

ur 

friction  velocity  'f  r  /  p 

r  * 

fl  -  (a/x)*  1 

8  = 

1 - 1  1 

*— • 

I 

P 

CO 

1—* 

o 

*r(a',/9) 

(a',/9) 

Bf  <a',/V) 

Br(a'./J) 

r 

f  («'-  1)  .(p-1)  , 

/a  (1  -o')  do 

J  0 

Dimensionless 

.  U 

qw 

=  heat  flux  parameter 

PCp  u fcji 

u+ 

dimensionless  velocity  in  boundary  layer  ,  defined  by  eq.  4.5 

x+ 

distance  along  the  wall,  defined  by  eq.  4.3 

xa,b 

distance  between  points  x  =  a  and  x  =  b,  defined  by  eq.  4.27 
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+ 

+ 


a 

e 


+ 


e 

e  . 

a,b,  c 

? 

Pr 

Sc 

Nu 


Re 


Re 


D 


St 


St 

Cf/2 


distance  normal  to  the  wall,  defined  by  eq.  4.4 

thermal  diffusivity,  defined  by  eq.  4.7 

eddy  viscosity,  defined  by  eq.  4.6 
T  -  T 

_ i_ 

T  -  T 
w  i 

value  of  6  at  point  C  for  constant  heat  flux  parameter  between 
x  =  a  and  x  =  b 

distance  from  the  wall,  defined  by  eq.  4.8 
Prandtl  number 
Schmidt  number 

hx 

Nusselt  number  — 

Reynolds  number  U’X 
v 

Reynolds  number  ^ 

Stanton  number  - ^-7; — 

pui  Cp 


Spalding  function 


friction  factor 


P  u 


— 2 

pu 


for  pipe  flow 


WATER  TANK 
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